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ABSTRACT 

We have used maximum entropy eclipse mapping to recover images of the visual 
surface brightness distribution of the primary component of the RS CVn eclipsing 
binary SV Cam, using high-precision photometry data obtained during three primary 
eclipses with STIS aboard the Hubble Space Telescope. These were augmented by 
contemporaneous ground-based photometry secured around the rest of the orbit. The 
goal of these observations was to determine the filling factor and size distribution of 
starspots too small to be resolved by Doppler imaging. The information content of 
the final image and the fit to the data were optimised with respect to various system 
parameters using the % 2 landscape method, using an eclipse mapping code that solves 
for large-scale spot coverage. It is only with the unprecedented photometric precision of 
the HST data (0.00015 mag) that it is possible to see strong discontinuities at the four 
contact points in the residuals of the fit to the lightcurve. These features can only be 
removed from the residual lightcurve by the reduction of the photospheric temperature, 
to synthesise high unresolvable spot coverage, and the inclusion of a polar spot. We 
show that this spottedness of the stellar surface can have a significant impact on the 
determination of the stellar binary parameters and the fit to the lightcurve by reducing 
the secondary radius from 0.794±0.009R Q to 0.727±0.009R Q . This new technique 
can also be applied to other binary systems with high precision spectrophotometric 
observations. 
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1 INTRODUCTION 

Star-spots and other forms of magnetic activity are preva- 
lent on rapidly rotating stars with temperatures low enough 
to have outer convective zones. The dynamo activity in these 
rapidly rotating cool stars leads to the suppression of con- 
vection over large areas of the stellar photosphere. These 
giant starspots can modulate the light of such active stars 
as they rotate by up to tens of percent. In general, stars 
that exhibit this type of activity are solar-type stars with 
rotational periods of less than one day. This rapid rotation 
is seen in half the G and K dwarfs in open clusters younger 
than 100 Myr, and may also persist into middle age through 
tidal locking in a close binary systems such as RS CVns. 

Doppler imaging is a powerful tool for determining 
where spots congregate on the stellar surface, but it is less 
successful at determining how much overall spot activity 
is present. Other methods such as TiO-band monitoring 



studies JO'Neal et al.lll99Sl) indicate that between 30-50% 
of the stellar surf ace is spotted comp ared to 10-15% from 
Doppler imaging (Jeffe rs et alJ 12002) . The discrepancy in 
these methods can be accounted for if the stellar surface 
is heavily covered by spots that are too small to be re- 
solved through Doppler imaging. To try to resolve this we 
have used the photometric capabilities of the HST (signal- 
to-noise ~ 5000) to eclipse-map the inner face of the F9V 
primary of the eclipsing bin ary SV Cam. Using these obser- 
vations, Ijeffers et alJ i2005l) found that the surface flux in 
the low-latitude region to be approximately 30% lower than 
the best fitting phoenix model atmosphere. A possible cause 
of this flux deficit is that the primary star's surface is pep- 
pered with unresolvable spots. They also found that when 
the 30% spottedness is extended to the entire stellar surface 
there still remains a flux deficit, which could be explained 
by the presence of a large polar cap. 
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Figure 1. The comprising sections of the primary eclipse from 
each of the 3 HST visits, with an off-set of 5xl0~ lx ergs/s/cm 2 
for clarity. 



2.2 Ground Based Photometry 

The ground based photometric observations were obtained 
with the 0.93-m James Gregory Telescope (JGT) at the Uni- 
versity of St Andrews, using a Wright Instruments CCD 
camera mounted at the Cassegrain focus. The observations 
were made with a narrow-band filter whose 85 A (FWHM) 
bandpass is centred at 5300 A, close to the V band. Using an 
exposure time of 60 s we obtained a photometric precision of 
±0.006 mag in good photometric conditions. The overhead 
for the chip readout between 60 s observations is typically 
15 s, giving an overall sampling interval of 75 s. Observations 
were made over a two week period centred at the time of the 
HST observations. The weather permitted the observation 
of four full primary eclipses and three full secondary eclipses 
over this period. 



In this paper we use t he Maximum Entropy 
eclipse-mapping code DpTS llCollier Cameron) Il997t 
ICollier Cameron fc Hildltchl Il997fr to firstly solve for the 
system parameters of the lightcurve, and then reconstruct 
the surface brightness distribution of SV Cam's primary 
star. We also show how a synthetic polar cap and the star 
being peppered with spots below the resolution limits of 
eclipse mapping are essential parameters for obtaining a 
good lightcurve fit. 



2 OBSERVATIONS 

SV Cam (F9V + K4V, P=0.59d, vsini=102 kms" 1 , i=90°) 
is a totally ecli psing, s h ort-p eriod binary system originally 
identified by iGuthnickl |l929). Contemporaneous HST and 
ground based observations were obtained to eclipse map SV 
Cam. The ground based photometric observations were ob- 
tained to complete the light curve outside primary eclipse 
making it possible to determine any global asymmetry in 
the lightcurve due to the presence of any uneclipsed spots. 



2.1 HST Observations 

We observed three primary eclipses of SV Cam during three 
HST visits at intervals of 2 days between 2001 November 
1 to 5. A total of 9 spacecraft orbits were devoted to the 
observations. The G430L grating of Space Telescope Imag- 
ing Spectrograph (STIS) was used to disperse starlight over 
2048 CCD elements. The observations are summarised in 
Tableland the lightcurve for each visit is shown in Fig. 
The observations cover a wavelength range from 2900 A 
to 5700 A, with an exposure time of 30 s and a cadence of 
40 s. The exposure time was chosen as a trade-off between 
maximising signal-to-noise per exposure and minimising the 
time interval between observations. The total count over the 
entire detector is 9.4 xlO 7 electrons per exposure, with the 
photometric precision calculated to be 0.00015 magnitudes 
(S:N = 5000). The main source of systematic error, due to 
the subtle motion of the satellite, was removed by cross- 
correlating each frame with respect to a template spectrum, 
and then adding or subtracting a scaled number of counts. 
After correcting for this 'jitter' the RMS scatter was consis- 
tent with the error expected from photon noise. 



3 DATA PROCESSING 
3.1 Data Reduction 

Reduction of the HST data comprised calibrating the two 
dimensional CCD images and extracting one dimensional 
spectra using the standard STIS data reduction pipeline. 
A detailed description of this procedure is contained in the 
STIS handbook. This pipeline routine produces bias and 
dark subtracted, flat-fielded wavelength calibrated images. 
The auto-wavecal was disabled for these images, but instead 
two 'Guest Observer' wavecals were taken per orbit. Cos- 
mic rays were removed using an algorithm that identifies 
and rejects cosmic rays and other non repeatable defects by 
comparing successive frames. 

The ground-based photometric data were reduced us- 
ing jgtphot, a software package develope d for use with th e 
James Gregory Telescope at St Andrews feell et al.lll993T) . 
The resulting light curves are calculated as differential mag- 
nitude values with respect to the marked comparison and 
check stars. 



3.2 Time Correction of HST data set 

As the time of the extracted one-dimensional spectra that 
is output from the standard STIS data reduction pipeline 
(calstis) is in UT, it is necessary to apply a further time 
correction to this value. The heliocentric correction due to 
the motion of the Earth around the Sun, and due to the 
orbital motion of the satellite was determined through use 
of the 'odelaytime' routine. To convert UT into atomic time 
32 leap seconds were added, and to convert atomic time 
into terrestrial time 32.184 seconds were added, resulting 
in the correct terrestrial time value. Absolute precision in 
HST time is required as it is necessary to have the primary 
eclipse of the HST and of the ground-based observations 
exactly aligned. 

3.3 Interpolation of Data Sets 

In order to analyse the absolute spectrophotometry from 
STIS together with the ground-based differential photome- 
try, it was necessary to place the magnitudes derived from 
the two instruments on a common magnitude scale. We 
constructed a synthetic filter, Vhst, with a bandpass and 
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Visit Obs. Date Frame No UT Start UT End Julian Date (+ 52210) Exposure Time (s) Phase Range 

1 01 Nov 01 1-48 21:04:33 21:47:44 4.879394-4.908758 30 0.901-0.950 

49-107 22:24:18 23:17:41 4.909383-4.971851 30 0.952-1.058 

108-165 00:00:23 00:51:58 5.001505-5.037328 30 1.107-1.168 

2 03 Nov 01 166-213 14:43:36 15:26:47 6.614915-6.644280 30 0.828-0.877 

214-272 16:02:23 16:55:46 6.669628-6.706701 30 0.920-0.983 

273-330 17:38:27 18:30:02 6.736343-6.772167 30 1.033-1.093 

3 05 Nov 01 331-378 14:43:36 15:26:47 8.416575-8.446561 30 0.865-0.916 

379-437 16:02:23 16:55:46 8.470752-8.507825 30 0.957-1.019 

438-495 17:38:27 18:30:02 8.537468-8.573291 30 1.069-1.13 



Table 1. Journal of observations for SV Cam using the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope 
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Figure 2. Interpolated HST data (primary eclipse) set plotted 
with JGT data set (secondary eclipse). 



4 LIGHTCURVE FITTING 

The combined HST and JGT lightcurve of SV Cam is fitted 
using the maximum entropy code DOTS. Included in DOTS is 
the ability to incorporate the surface geometry and radial ve- 
locity variations of tida lly distorted close binary components 
llCollier Cameronlll997l) while solving for starspot coverage. 

4.1 Eclipse Mapping 

By using eclipse mapping it is possible to determine the de- 
tailed locations of low latitude spots on the stellar surface. 
If spots are present on the inner surface of the primary, that 
is occulted by the cooler secondary star, jagged discontinu- 
ities will be produced on the eclipse profile as the primary is 
occulted by the secondary. The timescales and amplitudes 
of these discontinuities will reflect the distribution of spot 
sizes on the inner surface of the primary. 



effective wavelength matching those of the filter used for 
the JGT observations, (bandpass of 80A and effective wave- 
length 5300A) and applied it to the HST observations. Simi- 
larly we constructed a digital version of the Johnson B filter, 
Bjj S t. 

At the time of each JGT flux measurement during pri- 
mary eclipse we interpolated the Vhst magnitude and the 
{Bhst — Vhst) colour index from the two nearest points 
bracketing the same phase. The colour equation relating the 
two magnitude systems was then derived by a linear least- 
squares fit to the plot of SVj g t — Vhst versus {Bhst ~ Vhst), 
as shown in Fig. [5] This yielded the colour equation 

5V jgt - V h st = a(B h st - Vhst) + P, (1) 

where the coefficient of the colour term a — —0.0725392 ± 
0.05836 and the zero-point (5 = -8.69367 ± 0.03318 magni- 
tudes. 

Using this linear calibration, all of the Vhst magnitudes 
were converted to the JGT system. This procedure takes ac- 
count of the zero point off-set, the effect of the Vhst and the 
JGT observations having slightly different effective wave- 
lengths, and the different spectral responses of the two filter 
sets and the two CCDs. The interpolated HST data set is 
shown in Fig. |5|along with the JGT data. 



4-1.1 Maximum Entropy 

The maximum entropy method allows us to determine the 
simplest image of the primary's surface (in terms of its infor- 
mation content as quantified by the Shannon- Jaynes image 
entropy), that can reproduce, at a specified goodness of fit, 
the perturbations to the primary eclipse. 

For images where the mapping parameter fi is bolomet- 
ric surface brightness and so can take any positive value, the 
entropy takes the form, 

S = VC/i-mi-Ziln— ) (2) 

mi 

i 

Here m; is the default value that a pixel will have when 
there are no other constraints imposed by the data. In this 
analysis a restricted form of the entropy is used as the filling- 
factor model has been adopted. This combines the entropy 
of the spot image fi and of the photospheric image (1 — fi): 

S = ^(-/,lnA_(i_ /i )l ri il^)). (3 ) 
^-^ rrii (1 — rm) 

i 

To construct the final image, the values for the spot and 
photospheric fi are iteratively adjusted to maximise; 

Q = S - A X 2 (4) 

This is equivalent to maximising the entropy 5* over the 
surface of a hyper-ellipsoid, of constant \ 2 , m image space. 
This is bounded by the constraint surface at some fixed value 
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Table 2. Summary of observed stellar parameters for SV Cam from the literature where T pr i and T sec are the primary and secondary 
photospheric temperatures (K), R pr i and R se c are the primary and secondary radii (R©), i is the binary system inclination (degrees), q 
is the secondary:primary mass ratio and S/P indicates spectroscopic or photometric observations. 



of \ 2 - The Lagrange multiplier, A, is set so that the final so- 
lution lies on a surface with \ 2 — M, M being the number 
of measurements in the data set. By setting a Lagrangian 
multiplier, it is possible to determine the extrema (maxi- 
mum) of the entropy, the image with the least information 
content, subject to the constraint of obtaining the best fit 
to the data. A more detailed explanation of this, and other 
methods, is discussed by Collier Cameron (2001). 

It is important not to over fit x 2 during the reconstruc- 
tion, as this can result in an ill-fitting photometric fit and 
artifacts on the final spot map. The number of iterations was 
set to 25, as through rigorous testing we found this to be the 
value where the final image did not contain any distortions 
due to ill fitting. 

4.2 Geometric Parameters 

The success of a surface image reconstruction depends 
largely on a correct determination of the geometric parame- 
ters of the binary system (Vincent et al 1993). Independent 
reconstructions for the starspot parameters and the geomet- 
ric parameters of the system, are susceptible to the distor- 
tions in the lightcurve yielding the wrong system parameters 
or vice versa. If the system parameters are wrong, in general 
a satisfactory fit to the data is only obtained by increasing 
the amount of structure in the stellar surface brightness dis- 
tribution, usually leading to a greater total spot area than 
is present when the correct geometric parameters are used. 

The principal parameters that affect the shape of the 
light-curve of an unspotted star are the relative tempera- 
tures, the radii of the two stars and the system inclination. 
The relative surface brightnesses of the two components de- 
termine the relative depths of the primary and secondary 
eclipses. The radii of the two stars and the binary system 
inclination are related to each other by the orbital phases at 
which the four contact points occur (phase 0.906,0.977,1.023 
and 1.093 respectively). These relations are further compli- 
cated by the presence of starspots which in general have the 
effect of altering the depth of the eclipses and making the 
contact points asymmetrical. 

SV Cam has been the subject of numerous photometric 
and spectroscopic studies. Despite this, there still does not 
yet exist a established and reliable set of stellar parameters. 



The various published combinations of stellar parameters 
that can be found for SV Cam in the literature are sum- 
marised in Table [5] A discussion of the individual merits of 
each set of binary system parameters is beyond the scope of 
this paper. Given the wide range of geometric parameters 
available for SV Cam, it was necessary to determine a set of 
parameters that would give an optimal fit to our lightcurve. 

In fitting the photometric lightcurve of SV Cam we use 
the mass ratio (q=0.641±0.007), semi-amplitude velocity 
(Ki=121.86±0.76, K 2 = 190.17±1.73) and sys tematic radial 
velocity (-9.31±0.78) of lRucinski et al.l (120021) as their solu- 
tion for SV Cam's spectroscopic orbit is very well defined. 

4.3 Temperature 

Ijeffers et al.l i2005l) determined the temperature of the sep- 
arated prima ry and secondary components of SV Cam us- 
ing phoenix iAllard et al.l200Cll) model atmosphere spectra. 
Model atmospheres, ranging in temperature from 5600 K to 
6500 K in 100 K steps, were fitted to the spectrum of the pri- 
mary and secondary stars using \ 2 minimisation. The min- 
imum x 2 value corresponds to a temperature of 6013±19K 
and 4804±143K for the primary and secondary stars re- 
spectively. The primary star was then isolated and fitted 
using the same method, where the minimum temperature 
is 6038±58 K. In this work we use the independently fitted 
primary temperature as the error estimate is more reliable. 

4.4 Radii 

The radii for the primary and secondary stars were deter- 
mined using the robust grid search method. This method 
uses a grid of radii as input to DoTS ranging from 1.17 R© 
to 1.26 Kq for the primary star, and 0.73 Kq to 0.82 R© for 
the secondary star in 0.005 R© intervals. For each grid point 
a model with the specified set of geometric parameters was 
iterated 25 times to obtain the lowest value of the reduced 
X 2 . The inclination was not included as it scales as (sin i) — 1 , 
implying a negligible difference between inclinations of 90° 
and 85°. 

The results for the x 2 minimisation are shown in Fig. |3] 
in the form of a contour plot. The minimum x 2 value occurs 
at 1.238 ±0.007 R Q and 0.794 ±0.009 R Q for the primary 
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Figure 3. Contour plot of the \ 2 landscape for the primary and 
secondary radii. From the centre of the plot, the first contour 
ellipse is the 1 parameter lcr confidence limit at 63.8%, the second 
contour ellipse is the 2 parameter lcr confidence limit at 63.8%, 
and the third contour ellipse is the 2 parameter 2.6<r confidence 
limit at 99%. 
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Figure 6. The final spot map of SV Cam (HST + JGT data). 
Note that phase runs in reverse to longitude. Spots reconstructed 
us ing y 2 and loc ated at the quadrature points have been shown 
bv lJeffersl 120051) to be spurious. Note that phase runs in reverse 
to longitude. 
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Figure 5. Observed minus Computed residuals with no polar cap 
(HST data only). The four contact points are shown as vertical 
lines. 

and secondary stars respectively. The best fitting stellar pa- 
rameters for the combined JGT and HST data sets are sum- 
marised in Table 

4.5 Lightcurve Fit and Surface Brightness Image 

The resulting Maximum Entropy fit to the lightcurve, us- 
ing the stellar parameters as previously defined is shown in 
Fig. 3] The observed data minus computed lightcurve resid- 
uals are shown in Fig. [S] in enlarged form covering just the 
HST primary eclipse data. The stellar surface image that 
results from these fits is shown in Fig. |S] 

The model lightcurve of the best fitting binary system 
parameters is shown in Fig. 2] When the modelled light- 
curve is subtracted from the observed data the residuals are 
symmetric about phase 1, and have two large discontinuities 
that occur at the first and fourth contact points either side 
of the primary eclipse (Fig. |SJ . Given the location of these 
discontinuities, an obvious explanation would be that the 
radius of either of the two stars has been incorrectly deter- 



Paramctcr sol. 1 sol. 2 



Ephemeris 52214.34475 (MJD) 

T pri 6038 ±58 K 5935 ±38 K 

T 3 „ 4804 ±143 K 4804 ±143 K 

R pri 1.238 ±0.006 R 1.235 ±0.007 R 

H, ec 0.794 ±0.006 R 0.727 ±0.006 R 

Polar Cap 46.7° 



Table 3. Geometric binary system parameters computed for SV 
Cam. For s olution 1 , the prima ry and secondary temperatures are 
taken from Jcffe rs et alJ 120051) . while solution 2 shows the radii 
and polar cap size solved using a reduced photosphcric tempera- 
ture to synthesise the star being peppered with small spots. 

mined. However, the accuracy of our determination of the 
stellar parameters, as described in Section 4, can eliminate 
this possibility. 

The surface brightness distribution (Fig. |SJ| shows two 
large spot features at the quadrature points similar to 'active 
longitudes' that hav e been observed using photometry on 
oth er RS CVns (e.g. lOlah et alJ <199lft . lLanza et alJ (1200 ll) 
and lLanza et alJ i2002l) ). These spots are reconstructed at 
the equator as there is no eclipse information to determine 
the latitude. The spot at 100° is stronger as that quadrature 
of the star has a higher spot coverage, which is also evident 
in the lower light level of the photometric lightcurve at phase 
~0.75. The spots reconstructed at 25°, 335°, and 360° result 
from spots on the primary star that have been occulted by 
the secondary star during primary eclipse. The symmetry is 
due to the inability of DoTS to determine the hemisphere 
of the spot feature. 



5 UNRESOLVABLE SPOT COVERAGE 

5.1 Reduced Photospheric Temperature 

Ijeffersl i2005T) tested the limitations of surface bright- 
ness images from photometric data by modelling many 
sub-resolution spots on an immaculate SV Cam. Surface 
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Figure 4. Combined HST and JGT lightcurves, with a Maximum Entropy fit. 



brightness distributions reconstructed from these synthetic 
lightcurves show distinctive spots at the quadrature points. 
The presence of similar spot features in Fig. |S| could in- 
dicate that the surface of SV Cam is peppered with many 
spots that are below the resolution capabilities of the eclipse 
mapping technique. This would then be consistent with the 
results of TiO band monitoring studies IIO'Neal et al ] |!998ft 
which show that between 30% and 50% of a star's surface 
may be spotted. 

Further evidence for the peppering of t he SV Cam's 
prima ry star by small starspots is shown by Ijeffers et alJ 
(2005). As previously described in section 4, they deter- 
mined the temperatures of the two stars using best-fitting 
phoenix model atmospheres. It was also found that the pri- 
mary's surface flux is approx 28% lower than predicted by a 
phoenix model atmosphere at the best fitting effective tem- 
perature. Even taking into account the spot distributions 
as shown in Fig. |S| this flux deficit can only be accounted 
for if the primary star's surface is peppered with unresolved 
spots. As these spots are not resolvable using the eclipse- 
mapping technique they can lead to a under estimation of 
the star's photospheric temperature and will have the ef- 
fect of decreasing the flux deficit during the eclipse at all 
wavelengths. 

To determine the reduced photospheric temperature, re- 
sulting from the presence of ma ny unresolvable spots, we ex- 
tend the work of Ijeffers! i2005T) . An extrapolated solar spot 
size distribution is applied to an immaculate SV Cam, for 
1.8%, 6.1%, 18%, 48% and 100% area filling factors of spots, 
with T ph =6038 K, and T SJ ,=4538 K. Each of these spot dis- 
tributions is modelled as a photometric lightcurve, and is 
then used as input to the Maximum Entropy eclipse map- 
ping code. To obtain a satisfactory fit to the lightcurve it 
was necessary to reduce the photospheric temperature to 
values as shown in Fig. |7| A quadratic fit to these points 
gives a temperature for 28% spottedness of 5935 K which is 
the average temperature over hotter and unresolvable cooler 
temperature regions. 
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Figure 7. The decrease of primary star's photospheric tempera- 
ture as a function of percentage of spot coverage on its surface. 



5.2 Polar Spot 

The work of I.Teffers et alJ l|2005h also showed that in addi- 
tion to the 28% spot coverage in the eclipsed region of the 
primary star, there is an additional 12.5% flux deficit in 
non-eclipsed regions of the star. The additional flux deficit 
can only be explained by the presence of a polar cap which 
would have the effect of increasing the depth of the primary 
eclipse, and would lead to an incorrect lightcurve solution. 

We include a polar cap in our lightcurve solution by ex- 
tending the 2-dimensional grid search described in section 
4 to a 3-dimensional grid search including the polar spot 
size from 35° to 50°. The polar spot is assumed to be at 
4500 K, circular, centred at the pole, and is in addition to the 
reduced photospheric temperature as described above. For 
each polar spot size the minimum primary and secondary 
radii are determined using a \ 2 contour map. These mini- 
mum x 2 values are plotted as a function of polar spot size in 
Fig- El The best fitting polar cap size, 46.7° was determined 
from the minimum of a quadratic function fitted to these 
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Figure 8. Quadratic fit to the variation of \ 2 as a function of 
polar spot size. 
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Figure 9. Contour plot of the x 2 landscape for the primary and 
secondary radii with a 46.69° polar spot. From the centre of the 
plot, the first contour ellipse is the 1 parameter la confidence 
limit at 63.8%, the second contour ellipse is the 2 parameter la 
confidence limit at 63.8%, and the third contour ellipse is the 2 
parameter 2.6c confidence limit at 99%. 
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Figure 11. The final spot map of SV Cam (HST + JGT data) 
with two 46.7° polar caps and reduced photospheric temperature. 
Note that phase runs in reverse to longitude. Spots reconstructed 
us ing y 2 and lo cated at the quadrature points have been shown 
bv ljeffers! l2005h to be spurious as they can be reconstructed from 
models of high starspot coverage. Note that phase runs in reverse 
to longitude. 
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Figure 12. Observed minus Computed residuals for the case with 
a reduced photospheric temperature and a polar cap of 46.69° 
radius. The four contact points are shown as vertical lines. 



points. This value is in good agreeme nt with the indepen - 
dently determined value of 42±6° bv Ijeffers et alJ |2005). 
The grid search of radii is then repeated using a fixed value 
for the polar spot size. The results for the x 2 minimisation 
are summarised in Table 3 and are shown as a contour map 
in Fig. EH 

5.3 Final Lightcurve Fit 

The best fitting stellar parameters are used to create a Max- 
imum Entropy lightcurve fit (Fig. I TUP and stellar surface 
brightness image for the primary star (Fig. 1110 . The ob- 
served data minus computed lightcurve residual for the pri- 
mary eclipse, as shown in Fig. 1121 is significantly flatter than 
the residuals shown in Fig. |K| The secondary eclipse shows 
also a better fit than in Fig^] The reduced \ 2 value after 25 
iterations was 2.05, which is slightly reduced from the value 
of 3.55 obtained without a polar spot and a reduced pho- 
tospheric temperature. This result shows that the presence 



of a polar cap and sub-resolution spots have a significant 
influence on the O-C residuals. 

The Maximum Entropy surface brightness distribution, 
as shown in Fig. 1111 is similar to that shown in Fig. |S| The 
spot feature at 0° and 360° is not an artifact as it is in the 
eclipse path of the secondary, but results from time variable 
structure in the base of the primary eclipse. It is a symmet- 
ric spot feature as image reconstructions from photometry 
are not able to determine a spot's latitude. The spot fea- 
ture at 25° and 335° in Fig. |S| has disappeared due to the 
correct determination of the primary star's temperature and 
consequently the secondary star's radius. 



6 DISCUSSION AND CONCLUSIONS 

We have used eclipse mapping, based on the maximum en- 
tropy method to recover images of the visual surface bright- 
ness distribution of the primary component of the RS CVn 
eclipsing binary SV Cam. It is only with the unprecedented 
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Figure 10. Combined HST and JGT lightcurves, with a Maximum Entropy fit with a polar cap and reduced photospheric temperature 
included. 



photometric precision of the HST data that it is possible to 
see strong discontinuities at the four contact points in the 
residuals of the fit to the lightcurve. These features can only 
be removed from the O-C lightcurve by the reduction of the 
photospheric temperature, to synthesise high unresolvable 
spot coverage, and the inclusion of a polar spot. 

In fitting SV Cam's lightcurve we used the indepen- 
dently determi ned primary and sec ondary temperatures, as 
determined by Ijeffers et alj 1120051) . as input to our solu- 
tion. The primary and secondary radii were then deter- 
mined using a 'grid search' method. Although the 'grid 
search' method is computationally intensive, it does pro- 
vide reliable results. In contrast to this is the 'downhill sim- 
plex' amoeba method jPress et alJll992T) which repeatedly 
returned lightcurve solutions from local \ 2 minima rather 
than from the global \ 2 minimum. The resulting parame- 
ters (R pr ;=1.238R0 an d R se c=0.794R(T)) are in closest agree- 
ment with the values of lPoimanskil ll998l) derived from spec- 
troscopic data (R pr i=1.25R© and R S ec=O.8R0). The errors 
derived in this work are an order of magnitude smaller than 
previous lightcurve solutions. 

The Maximum Entropy fit to the data using the best- 
fitting orbital parameters shows strong discontinuities at the 
four contact points of the primary eclipse. When we recon- 
structed an image using only the ground based JGT data 
(photometric precision = 0.006) we found that these strong 
discontinuities were not visible in the residual lightcurve. 
It is only with the high photometric precision of the HST 
(0.00015 mag. or S/N 5000) that it is possible to distinguish 
the strong discontinuities in the residual lightcurve. 

The surface brightness distribution shows spots at the 
quadrature points consistent with surface images of other RS 
CVns recon structed using y 2 minimisation e . g. lOlah et alJ 
ll99ll) and lLanza et all 1200 j) . lLanza et all 12002?) . Ho~ 
ever, the relia bility of these images has been investigated by 
Ijeffersl l2005h where spurious spots at the quadrature points 
have been reconstructed from a synthetic star containing a 



high degree of sub-resolution spots. This is consistent with 
results from other methods such as TiO-band monitoring 
lO'Neal et alJll99d) which indicate that between 30% and 
50% of a star's surface could be covered in starspots at all 
times. 

A high total spot coverage on SV Cam's primary star 
will modify t he apparent photosp heric temperature. In a 
related paper Ijeffers et alJ l2005h showed that when the 
surface flux in the low-latitude eclipsed region was ap- 
proximately 30% lower than the the best fitting phoenix 
model atmosphere. In this paper we showed that this is 
equivalent to a reduction in the photospheric temperature 
from 6039±58 K to 5935±38 K. High spot coverage also 
has important s tructu ral implications. As investigated by 
ISpruit fc Weissl ll98rJl . a star with high spot coverage will, 
over thermal timescales, readjust its structure to compen- 
sate in radius and temperature. For example the primary 
star's radius is 10% larger than expected for its spectral 
type, and could provide an explanation as to the large vari- 
ation in binary system parameters as shown in Table |5] 

The addition of a polar spot and the reduction of 
the photospheric temperature has a negligible influence on 
the empirically determined radius of the primary star, but 
a more significant influence on the secondary, decreasing 
it from 0.794±0.007R Q to 0.727±0.009 R . The presence 
of a polar spot will increase the depth of the primary 
eclipse, while a reduction in the photospheric temperature 
will decrease the depth of the eclipse. In both cases the 
binary system compensates for this by making the sec- 
ondary star smaller as it consequently needs to eclipse less 
light. The standard luminosity-radius-temperature relation 
(L=47rR 2 aT 4 ) shows that if T pr i and T sec are fixed the 
only variable parameters are the radii of the primary and 
secondary stars. The timing of the first and fourth contact 
points fixes the primary radius, leaving the secondary radius 
as the only adjustable parameter. 

The resulting residuals from the Maximum Entropy 
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lightcurve fit are significantly flatter with the lightcurve so- 
lution that includes a polar cap and a reduced photospheric 
temperature. However, there are still variations in the resid- 
ual lightcurve, which could result from the eclipse mapping 
routine interpreting small scale variations in the lightcurve 
(i.e. unfitted spots) as noise. A 0.005 change magnitude as 
shown in Fig. II 21 would result in additional 14% of the area 
of the annulus around the secondary being spotted than is 
shown in the surface brightness images. However, the re- 
maining structure in the O-C residual could result from an 
incorrect value of the limb darkening or the presence of plage 
on the primary star's surface, which will require further in- 
vestigation. 

The surface brightness image reconstructed with the re- 
fined binary system parameters (Fig. M H has reduced sur- 
face structure compared to the image reconstructed with 
the original parameters. This is noticeable at longitudes 25° 
and 335°, where the small spots have disappeared with the 
refined binary system parameters. The structure at longi- 
tudes 0° and 360° are not artifacts of the eclipse mapping 
technique as they are in the eclipse path of the secondary, 
but result from time variable structure on the surface of the 
primary star during primary eclipse. The features are sym- 
metric as it is not possible to resolve a starspots latitude 
using photometric observations. Doppler images of active 
stars show that polar caps do not have a uniform structure, 
unlike the polar caps used in our model and shown in Fie. 1111 

The upcoming COROT and KEPLER missions, which 
are designed to detect transits eclipses of stars by terres- 
trial sized planets, will discover thousands of eclipsing bi- 
nary stars with micromag photometry, while missions such 
as GAIA should deliver 10 6 new eclipsing binary stars. The 
benefit of such a wealth of high precision data will be lost 
if it is not possible to accurately solve the binary system 
parameters from the lightcurves. 
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